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Background 

The  recent  years  have  seen  numerous  examples  where  mechanics,  in  conjunction  with  systematic 
mathematical  analysis  and  large-scale  computation,  has  played  a  critical  role  in  the  discovery, 
development  and  improvement  of  new  materials.  I  he  discovery7  of  extra-ordinarily  large  strain 
in  ferromagnetic  shape-memory  alloys,  the  discovery  of  new  shape-memory7  alloys  with  ultra- 
low  hysteresis,  morphological  instabilities  in  thin  films  and  development  of  new  quantum  dot 
structures  are  but  a  few  examples.  However,  each  problem  has  created  its  own  small  community 
and  the  endeavors  have  been  largely  disconnected  from  each  other.  With  the  growing  list  of 
successes  and  the  growing  communities,  the  time  is  ripe  to  have  a  systematic  discussion  about 
creating  a  broader  approach  wherein  the  diverse  groups  can  learn  from  each  other  and  new7 
problems  can  be  tackled.  This  symposium  successfully  created  a  forum  for  this  discussion. 

Symposium 

This  symposium  highlighted  how  mechanics  across  scales  ranging  from  angstroms  to  meters  has 
enabled  the  improvement  of  existing  materials  and  the  discovery7  of  new  materials.  It  brought 
together  speakers  from  diverse  disciplines  including  mechanics,  mathematics  and  materials 
science  and  backgrounds  ranging  from  academics  to  industry7  and  working  on  diverse  problems. 
The  breadth  of  this  symposium  can  be  gauged  by  the  range  of  exciting  topics,  including: 

•  Texture  and  grain  growth 

•  Thin  film  mechanics 

•  Fundamental  issues  in  elasticity,  including  symmetry,  the  Eshelby  problem,  and 
minimum  principles  for  rate  problems 

•  Control  devices  using  shape-memory  alloys 

•  Covalent  bonding  in  biological  systems,  the  mechanics  of  viruses  and  human  organs,  and 
the  use  of  mechanics  in  vascular  NiTi  stent  design 

•  Homogenization  theory  for  dynamic  linear  continuum  systems  and  shape-memory 
honeycomb  structures 
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•  Experiments,  modeling  and  simulation  on  the  mechanics  of  advanced  materials,  including 
carbon  nanotubes,  shape-memory  alloys,  thermoplastic  elastomers,  ferromagnetic  shape- 
memory  alloys,  and  piezoelectrics 

•  Analysis  and  application  of  quasicontinuum  methods 

The  talks  described  various  examples  and  contemporary  methods  of  analysis,  computation  and 
experiments  that  have  been  brought  to  bear  on  this  subject,  and  set  the  stage  for  broad-ranging 
discussions. 

This  symposium  was  organized  as  a  part  of  the  45tn  Annual  Technical  Meeting  of  the  Society  of 
Engineering  Science  and  was  held  in  honor  of  Richard  D.  James,  winner  of  the  2008  Prager 
Medal. 

The  list  of  speakers  is  in  Appendix  A,  and  the  abstracts  that  have  been  extracted  from 
Conference  Materials  in  Appendix  R. 

Disbursement  of  Funds 

The  speakers  are  active  and  recognized  researchers  in  the  field,  and  therefore  the  organizers 
anticipated  that  a  majority  of  them  will  be  able  to  obtain  other  sources  of  funds  to  attend  this 
symposium.  At  the  same  time,  organizers  of  the  symposium  sought  to  bring  together  researchers 
at  various  stages  of  their  careers,  in  particular  to  draw  from  the  experience  of  very  senior 
researchers  and  to  inspire  a  selected  group  of  early-career  researchers.  The  AFOSR  funds  were 
used  to  support  the  hotel  and  registration  expenses  of  5  early-career  researchers,  and  to  support 
the  registration  expenses  of  1  senior  researcher. 
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Appendix  A:  List  of  Speakers 

Symposium  on  Using  Mechanics  to  Discover  New  Materials 

45th  Annual  Technical  Meeting  of  the  Society  of  Engineering  Science 
October  12-15,  Urbana-Champaign 


In  order  of  presentation: 

1.  David  Kinderlehrer,  Carnegie  Mellon  University 

2.  Irene  Fonseca,  Carnegie  Mellon  University 

3.  Scott  J.  Spector,  Southern  Illinois  University 

4.  Raffaella  Rizzoni,  University  of  Ferrara 

5.  Liping  Liu,  University  of  Houston 

6.  L.  Ben  Freund,  Brown  University 

7.  William  S.  Klug,  University  of  California  Los  Angeles 

8.  Isaac  Chenchiah,  University  of  Bristol 

9.  Rodney  Clifton,  Brown  University 

10.  John  Willis,  University  of  Cambridge 

1 1 .  Graeme  Milton,  Utah  University 

12.  Michael  Ortiz,  California  Institute  of  Technology 

13.  Traian  Dumitrica,  University  of  Minnesota 

14.  Kaushik  Dayal,  Carnegie  Mellon  University 

15.  Mitchell  Luskin,  University  of  Minnesota 

16.  Pedro  Ponte-Castaneda,  University  of  Pennsylvania 

17.  Valery  Levitas,  Iowa  State  University 

1 8.  Ellad  Tadmor,  University  of  Minnesota 

19.  Anja  Schlomerkemper,  Max-Planck  Institute 

20.  Ichiro  Takeuchi,  University  of  Maryland 

21.  Doron  Shilo,  Technion 

22.  Ryan  Elliott,  University  of  Minnesota 

23.  Pradeep  Sharma,  University  of  Houston 

24.  Nicolas  Triantafyllidis,  University  of  Michigan 

25.  Timothy  Healey,  Cornell  University 

26.  Dominique  Schryvers,  Antwerp  University 

27.  Brian  Berg,  Boston  Scientific  Inc. 
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Appendix  B:  Presentation  Abstracts  (extracted  from  Conference  Materials) 
Symposium  on  Using  Mechanics  to  Discover  New  Materials 

45th  Annual  Technical  Meeting  of  the  Society  of  Engineering  Science 
October  12-15,  llrbana-Champaign 
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Prager  Medal  Symposium  in  Honor 
of  Richard  James:  Using  Mechanics 
to  Discover  New  Materials 


ORGANIZERS: 

Kaushik  Bhattachar  ya,  Division  of  Engineering  and  Applied  Science,  California  Institute  of  Technology 
Kaushik  Dayai,  Carnegie  Mellon  University 

Antonio  DeSimone,  SISSA 
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New  Perspectives  on  Texture  Evolution 


David  Kinderlehrer 

Professor 

Mathematical  Sciences,  Carnegie  Mellon  University 

Shlomo  Taasan 

Professor 

Mathematical  Sciences,  Carnegie  Mellon  University 

Nearly  all  technologically  useful  materials  are  polycrystalline,  consisting  of  small  crystallites,  called  grains,  separated  by 
interfaces,  the  grain  boundaries.  The  energetics  and  arrangement  of  this  network,  itstejxture,  are  important  factors  in 
material  behavior,  and  constitute  the  basic  problem  of  micro  structure.  Experimentally,  the  inverse  relationships  between 
grain  boundary  populations  and  energy  have  been  observed  formany  years  We  introduce  the  grain  boundary  character 
distribution  (GBCD),  a  new  statistic,  and  offer  compellingevidence  for  its  strong  depen  den  ce  on  grain  boundary  energy 
In  addition  we  demonstrate  that  there  are  natural  laws  for  texture  and  the  GBCD  usinglargenetworksimulations.  This  is 
joint  w  ork  with  the  MIMP  group  at  Carnegie  Mellon  University 
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Equilibrium  Configurations  of  Epitaxially  Strained  Crystalline 
Films:  Existence  and  Regularity  Results 

Irene  Fonseca 

Department  of  Mathematical  Sciences,  Carnegie  Mellon  University 

Strained  epitaxial  films  grown  on  a  relat  ively  thick  substrate  are  considered  in  the  context  ofplane  linear  elasticity.  The 
total  free  energy  of  the  system  is  assumed  to  be  the  sum  of  the  energy  of  the  free  surface  of  the  film  and  the  strain  aiergy. 
Due  to  thelattice  mismatch  between  film  and  substrate,  flat  configurations  are  in  general  energetically  unfavorable  and  a 
corrugated  or  islanded  morphology  is  the  preferred  growth  mode  of  the  strained  film  New  regularityresults  forvolume 
constrained  local  min  imizers  of  the  total  free  aiergy  are  established,  leading  to  a  rigorous  proof  ofthe  zero  contact  angle 
condition  betweai  islands  and  wetting  layers. 


ABSTRACTS 
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On  the  Symmetry  of  Energy  Minimizing  Deformations  in 
Nonlinear  Elasticity:  Compressible  Materials 


Scott  J,  Spector 

Professor  of  Mathematics 

Department  of  Mathematics,  Southern  Illinois  University 

Jeyabal  Sivaloganathan 

Reader  in  Mathematics 

Department  of  Mathematical  Sciences,  University  of  Bath 

Consider  a  homogeneous,  isotropic,  hyperelastic  body  that  occupies  a  three-dimensional,  thick  spherical  shell,  S,  in  its 
reference  state  and  is  subject  to  radially  symmetric  displacement  boundary  conditions  on  its  inner  and  outer  boundaries. 
It  is  shown  [1,2]  that,  to  each  deformation  of  the  shell,  there  corresponds  a  radial  deformation  that  has  less  elastic  energy 
than  the  given  deformation,  whenever  the  stored- energy  function  is  polyconvex  and  grows  sufficiently  rapidly  at  infinity 

The  key  ingredient  is  a  new  radial-symmetrization  procedure  that  is  appropriate  for  problems  where  the  symmetrized 
mapping  must  be  one-to-one  in  order  to  prevent  interpenetration  of  matter  The  radial  symmetrization  of  an  orientation 
preserving  diffeo  morph  ism,  u.S — >S*,  between  shells  S  andS*  is  the  radial  deformation 

v(x)=(R(r)/r)x,  where  r=|x|, 

that  maps  each  sphere  Sr  of  radius  r>0,  centered  at  theorigin  into  another  such  sphere  that  encloses  the  same  volume 
as  u(Sr).  Since  the  volumes  enclosed  by  the  surfaces  u  (Sr)  andv(Sr)  are  equal,  the  classical  isoperimetric  inequality 
then  implies  that  the  surface  area  ofu(Sr)  is  biggerthan  the  surface  area  ofv(Sr).  The  equality  ofthe  enclosed  volumes 
together  with  this  reduction  in  surface  area  is  then  shown  to  give  rise  to  a  reduction  in  total  energy  for  many  ofthe  consti 
tutive  relations  used  in  nonlinear  elasticity 

These  results  are  also  extended  to  dasses  ofSobolev  deformations  and  applied  toprove  that  the  radially  symmetric  solu 
tions  to  these  boundary  value  problems  are  global  energy  minimizes  in  the  dassof  (possibly  non  symmetric)  deforma 
tions  of  a  three-dimensional  thick  spherical  shell. 

References 

[1]  J.  Sivaloganathan,  S.J.  Spector,  Http://wvmmath.siu  edu/spector/incompressible.pdf 

[2]  J.  Sivaloganathan,  SJ.  Spector,  Http  ://www.  math  siu  edu/spector/compressible.pdf 
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Shape  Control  of  Polymeric  Sheets  Using  SMA  Wires 

Dr.  Raffaella  Rizzoni 

Department  of  Engineering,  University  of  Ferrara 

Mattia  Merlin 

Department  of  Engineering,  University  of  Ferrara 

We  investigate,  both  from  an  expenmmtal  and  a  theoretical  point  of  view,  the  behavior  ofSMA  composites  fabricated  by 
bonding  thin  pre-deformed  SMA  wires  to  a  polymeric  matrix.  The  motivation  of  this  study  is  primarily  the  design  ofa 
structure  changing  shape  to  change  operating  characteristics.  Some  relevant  material  properties  are  experiment  ally  deter 
mined  A  large  number  oftests  are  carried  out  on  Ti  SO. Sat. %Ni  wires,  inorderto  describe  the  stress-temperature  and  the 
strain  temperature  relations.  In  particular,  the  transformation  temperature  ranges  ofthe  wires  are  obtained  via  differential 
scanning  calorimetry  (DSC)  tests  and  the  latent  heats  on  heating  and  cooling  are  calculated.  The  elastic  moduli  ofthe 
martensite  and  austmite  phases  are 

det  erm  ined  by  toisile  tests  at  low  and  high  temperatures.  Constant  stress  tests  and  constrained  recovery  tests  are  also 
performed,  in  order  to  determine  the  relation  between  stress  and  the  transformation  temperatures,  and  the  maximum 
recovery  stress,  respectively.  The  mechanical  properties  of  the  polymeric  matrix  are  also  studied.  To  simulate  the  behavior 
ofthe  composite  a  simple  model  is  presented,  which  assumes  classical  plate  theory  for  the  composite.  The  behavior 
ofthe  SMA  wires  is  describedbythetheorv  for  heterogeneous  thin  wires  developed  in  [2],  which,  using  ideas  ofthe 
constrained  theory  of  martensite  [l],  is  specialized  to  phase  transforming  materials  with  high  elastic  moduli  [3].  Simula 
tions  ofthecomposite  changeofshape  for  different  arrangements  ofthe  SMA  wires  aregiven.  Experimmtal  results  on  an 
epoxy  resin  composite  plate  with  SMA  embedded  wires  are  also  presented,  showing  that  an  appreciable  bending  can  be 
obtained  during  wires  activation. 

References 

[  l]  J.M.  Ball,  C.  Chu  and  R.D  James,  (1995),  Hysteresis  during  stress- induced  variant  rearrangement,  Jour,  de  Physique  IV,  supple¬ 
ment  to  Journal  de  Physique  III,  Vol.5, 

colloque  C8, 245-251. 

[2]  H.  Le  Dret  and  N.  Meunier,  (2005),  Modeling  heterogeneous  martensitic  wires,  Math  Models  and  Meth  AppLSci.,  15,3, 
375—406. 

[3]  R.  Rizzoni,  Modelmgelastic  sheets  with  bondedshape  memory  wires,  preprint. 
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Solutions  to  the  Periodic  Eshelby  Inclusion  Problem  in  Two 
Dimensions 


Liping  Liu 

Postdoc 

California  Institute  of  Technology 

We  consider  the  classic  Eshelby  inclusion  problem  on  a  finite  unit  cell  with  periodic  boundary  conditions.  We  give  solu 
tions  in  terms  of  Cauchy  type  integrals  to  this  problem  in  two  dimensions  for  general  shapes  ofin  elusion  and  nonuni 
form  but  divergence  free  eigenstress.  This  representation  formula  is  reminiscent  ofthe  familiar  Greens  representation 
formula.  Mor^  it  is  shown  that  a  Vigdeigauz  structure  does  not  have  the  Eshelby  uniformity  property  for  non  dilatational 
uniform  eigenstress  unless  it  degenerates  to  a  simplelaminate. 


106 


SOCIETY  OF  ENGINEERING  SCIENCE  ♦  45TH  ANNUAL  TECHNICAL  MEETING 


Modeling  the  Forced  Separation  of  a  Molecular  Bond 

LB.  Freund* 

Division  of  Engineering  Brown  University 
freund@brown.edu 

Interest  in  this  topic  derives  from  the  observed  behavior  ofthenon  covalent  molecular  bonds  that  account  for  adhesion 
of  eukaryotic  animal  cells  to  extracellular  matrix  [l].  For  purposes  of  modeling,  bonding  is  understood  to  be  the  confine 
ment  of  the  configuration  of  the  the  system  within  an  energy  well,  which  arises  as  a  principal  feature  of  an  overall  energy' 
landscape  [2].  One  consequence  of  application  of  a  force  on  the  bond  is  to  distort  this  landscape  over  time.  The  system 
functions  in  a  thermal  environment,  and  the  depth  ofthe  bond  well  is  assumed  throughout  to  be  significantly  greater 
than  the  thermal  energy  unit  kT.  The  bond  state  is  characterized  by  a  probability  distribution  over  the  energy  landscape 
and,  as  the  landscape  distorts,  its  evolution  is  governed  by  a  partial  differential  equation.  The  system  is  analyzed  in  order 
to  extract  information  on  the  most  probable  time  for  bond  separation,  the  maximum  force  that  can  be  sustained  by  the 
bond,  and  the  sensitivity  ofbond  separation  behavior  to  loading  rate,  for  example.  The  goal  of  experiments  in  which 
molecular  bonds  a  re  separated  under  force  isto  infer  magnitudes  bond  characteristics.  The  purpose  ofthis  discussion  is 
to  consider  the  relationships  between  measured  quantities  and  properties  ofthe  bond  itself  as  modulated  by  an  int  erven- 
ing  loading  apparatus  [3]. 

References 

[1]  E.  Evans  and  K.  Ritchie,  Dynamic  strength  of  molecular  adhesion  bonds,  Biophysical  Journal  72, 1541-1555, 1997. 

[2]  P.  Hanngi,P.Tallcner  and  M.  Morlcovec,  Reaction-rate  theory  Fifty  years  after  Kramers,  Reviews  ofModern  Physics  62,  251-342, 
1990. 

[3]  L.  B  Freund,  Observations  on  the  mechanics  ofa  molecular  bond  under  force,  in  (it  Statistical  Mechanics  ofCellular  Systems  and 
Processes},  ed.  by  M  H.  Zaman,  Cambridge  University  Press,  in  press,  200S. 
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Modeling  the  Mechanics  of  Viral  Capsids 

William  S.Klug 

Depart  front  of  Mechanical  and  Aerospace  Engineering,  University  of  California,  Los  Angeles 

Melissa  M.  Gibbons 

Graduate  Student 

Department  of  Mechanical  and  Aerospace  Engineering,  University  of  California  Los  Angeles 

Lin  Ma 

Graduate  Student 

Department  of  Mechanical  and  Aerospace  Engineering,  University  of  California,  Los  Angeles 

As  revealed  by  techniques  of  structural  biology  the  protein  shells  of  viruses  (capsids)  are  some  of  natures  most  beauti¬ 
ful  and  robust  examples  of  highly  symmetric  multiscale  self-assembled  structures.  The  ability  ofviral  capsids  to  respond 
structurally  and  mechanically  to  physical  and  chemical  stimuli  also  gives  them  tremendous  potential  as  components  for 
the  design  of  synthetic  materials  with  adaptive  properties,  as  has  been  demonstrated  recently  by  the  creation  ofvirus- 
based  nanowires.  A  series  of  recent  indentation  experiments  using  atomic  force  microscopy  (AFM)  has  shown  that 
capsidsalso  possess  impressive  mechanical  properties  of  strength  and  elasticity.  In  this  talk  I  will  present  some  theory 
and  numerical  modeling  ofviral  capsids  based  on  continuum  mechanics  and  Ginzbuig-Landau  theory  ofphase  transi¬ 
tions,  and  I  will  discuss  what  weve  learned  about  the  elastic  response  and  mechanical  failure  ofviral  capsids,  and  coupling 
otglobal  capsid  mechanical  responseto  local  protein  conformational  change.  Overall  the  talk  will  highlight  ourongoing 
efforts  to  push  the  limits  of  usefulness  of  continuum  theoryvia  coupled  continuum  atomistic  multiscale  modeling  of 
capsids  and  other  laig^  protein  assemblies. 
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A  Model  for  T4  Tail  Sheaths 


Isaac  Chenchiah 

Lecturer 

Department  of  Mathematics,  University  of  Bristol 

Valeriy  Slastikov 

Lecturer 

Department  of  Mathematics,  University  of  Bristol 

Building  on  the  insights  of  Falk  and  James  [l]  we  present  a  model  for  bacteriophage  T4  tail  sheaths.  The  model  predicts 
helical  configurations,  phase  transitions  and  non  uniform  states — these  are  analysed  and  compared  with  experiments. 

Reference 

[l]  Falk,  Wand  James,  R.  D,,  Elasticity  theory  for  self-assembled  protein  lattices  with  application  to  the  martensitic  phase  transition  in 
bacteriophage  T4  tail  sheath,  Physical  Review  E  73(1),  011917  (2006) 
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Effective  Relations  for  Inhomogeneous  Elastodynamics  and 
Electrodynamics 
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Professor 
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J  R.Willis^da mtp.cam.ac.uk 

Thetheoryun  dedying  effective  const  it  ut  Re  relations  for  inhomogeneous  elastodynamics  is  largely  complete  even 
though  there  may  remain  practical  difficulties  in  implementation.  In  contrast,  there  appears  still  to  be  some  resort  to  ad 
hoc  prescription  in  the  case  of  inhomogeneous  electrodynamics,  even  in  the'homogenization  limit"  The  purpose  of  this 
talk  will  be  to  follow,  as  faithfully  as  possible,  a  route  recently  developed  and  explored  for  elastodynamics  [1,2],  in  the 
context  of  electro  dynamics,  to  establish  correspondences  and  differences  and  perhaps  improve  thereby  appreciation  of 
the  structure  of  electro  dynamic  effect  Re  relations. 
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Minimization  Variational  Principles  for  Acoustics, 
Elastodynamics,  and  Electromagnetism  in  Lossy 
Inhomogeneous  Bodies  at  Fixed  Frequency 
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The  classical  energy  minimization  principles  ofPirichlet  and  Thompson  are  extended  as  minimization  principles  to 
acoustics,  el asto dynamics  and  electromagnetism  in  lossy  inhomogeneous  bodies  at  fixed  frequency  This  is  done  by 
building  upon  ideas  ofCherkaev  and  Gibiansky,  who  derived  minimizat  ion  variational  principles  for  quasistatics.  In 
the  absence  of  body  forces  the  primary  el  asto  dynamic  minimization  variational  principles  have  a  minimum  which  is 
the  time -averaged  elastic  energy  dissipated  to  heat  in  the  body  The  variational  principles  provide  constraints  on  the 
boundary  values  of  the  fields  when  the  moduli  are  known.  Conversely  when  the  boundary  values  of  the  fields  have  been 
measured,  then  they  provide  information  about  the  values  of  the  moduli  within  the  body.  This  should  have  application 
to  tomography.  We  also  derive  saddle  point  variational  principles  which  correspond  to  variational  principles  ofGurtin, 
Willis,  andBorcea. 
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Minimum  Principles  for  the  Trajectories  of  Systems  Governed 
by  Rate  Problems 
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Theevolution  of  many  physical  systems  is  the  result  of  a  balance  bet ween  dissipative  and  energetic  forces.  Thisbalance 
results  in  evol  ut  i  on  a  ry  equations  that  govern  the  time-dependence  of  the  state  ofthe  system.  The  classical  treatment  of 
these  equations  tends  to  break  down  when  applied  to  strongly  nonlinear  systems  that  develop  evolving  microstructures. 
In  these  cases,  the  energy  ofthe  system  lacks  lower  semicontinuity  and  equilibrium  solutions  that  minimize  the  energy 
do  not  exist  in  general  In  particular,  time  incremental  formulations  suffer  from  a  restart  problem,  i.e.,  the  incremental 
problem  becomes  ill  defined  when  a  microstructure  has  been  previously  established  in  thesystem  and  the  initial  condi 
tions  for  the  future  evolution  ofthe  system  are  only  defined  in  the  sense  of  Young  measures.  However,  it  is  possible 
characterize  entire  trajectories  ofthe  system  as  minimizers  of  certain  energy  dissipation  functionals,  thus  entirely  cir 
cumventingthe  restart  problem  ofthe  classical  rate  and  time-incremental  problems.  Because  ofthe  minimizingproperty 
ofthe  trajectories  ofthe  system,  tools  ofthe  calculus  ofvariations  such  as  the  direct  method  for  establishing  existence, 
relaxation  and  opt  imal  scaling  can  be  applied.  Selected  examples  illustrate  how  those  tools  shed  light  into  microstructural 
evolution  and  on  the  effective  kinetics  governing  the  macroscopicevolution  ofthe  system. 
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Nanomechanics  of  Carbon  Nanotubes 


Traian  Dumitrica 

Professor 

Department  of  Mechanical  Engineering,  University  of  Minnesota 

Advances  in  condensed  matter  physics  and  theoretical  chemistry  have  made  possible  a  comprehensive  modeling  of  mate 
rials.  We  can  adapt  and  apply  these  theoretical  methods  to  study  the  am  azing  properties  of  nano  structures.  I  will  discuss 
two  examples  where  such  theories  helped  the  understanding  of  the  mechanical  response  of  carbon  nanotubes  (generally 
viewed  as  a  paradigm  for  nanoscale  materials). 

In  materials  modeling  the  most  common  wayto  implement  molecular  dynamics  is  via  translational  periodic  boundary 
conditions.  This  is  not  the  natural  choice  when  modeling  carbon  nanotubes.  Using  the  helical  and  rotational  symme¬ 
tries  of  the  nanoscale  graphitic  tubules,  molecular  dynamics  and  structural  relaxation  can  be  done  in  a  simplified  way, 
on  a  modest  number  of  atoms.  This  new  method  [l],  termed  objective  molecular  dynamics  [2,  3],  is  compatible  with 
quantum  mechanics  under  the  Bom  Oppenheimer  approximation.  Theutility  of  objective  molecular  dynamics  will  be 
presented  in  the  context  of  studying  the  carbon  nanotubes  under  bending  and  twist. 

Combining  a  probabilistic  approach  ofthe  rate  theory  with  detailed  quantum  mechanical  computations  of  failure 
nucleation  and  trail  sit  ion -state  barriers,  allows  for  a  comprehensive  analysis  ofthe  underlying  atomic  mechanisms  and 
evaluation  ofthe  yield  strain  for  arbitrary  nanotubes  under  realistic  conditions  [4],  The  numerical  results  are  captured  in 
a  concise  set  ofequations  for  the  breaking  strain,  and  reveal  a  competition  between  two  alternative  routes  ofbrittlebond 
breaking  and  plastic  relaxation.  The  employed  probabilistic  approach  ultimately  allows  for  the  creation  of  a  “strength 
map”,  which  plots  the  likelihood  that  a  nanotube  will  break — and  how  its  likely  to  break. 
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Phonon  Analysis  of  Carbon  Nanotubes  with  Arbitary 
Chirality  by  the  Objective  Structures  Framework 
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Department  of  Aerospace  Engineering  and  Mechanics,  University  of  Minnesota 

Normal  modes  of  vibrations  of  collections  of  atoms,  called  phonons,  provide  important  information  about  thermody¬ 
namic  properties  of  materials.  Further,  they  canbe  exploited  to  predict  equilibrium  crystal  structures  as  a  function  of 
temperature  and  other  environmental  variables  using  bifurcation  techniques  [1,2].  However,  computing  the  normal 
mode  spectrum  is  extremely  expensive,  except  in  the  special  case  of  periodic  crystals  where  Fourier  techniques  simplify 
the  calculation  Recently  R.  D.  James  observed  a  close  analogybetween  crystals  and  certain  complex  atomic  structures 
as  well  as  biological  structures  [l].  Structures  that  satisfy  this  analogy  are  termed  “Objective  Structures*  They  include 
carbon  nanotubes  ofall  chiralities,  buckyballs,  as  well  as  many  biological  virus  components.  This  analogy  brings  in  the 
idea  of“Objective  Bound  ary  Conditions*  The  Objective  Boundary  Conditions  allow  for  the  same  efficiencies  that  have 
made  crystal  line  structures  accessible  to  efficient  calculation  with  periodic  boundary  conditions  (ex.  [4]),  by  allowing 
treatment  of  a  single  unit  cell  representative  ofthe  entire  structure.  Here,  we  develop  the  Objective  analog  ofthe  Fourier 
transform  to  analyze  phonon  modes  in  carbon  nanotubes.  It  allows  us  to  treat  carbon  nanotubes  with  great  efficiency  by 
providing  a  basis  that  decouples  the  modes.  Further,  previous  calculations  of  phonon  spectra  apply  conventional  periodic 
boundary  conditions  along  the  length  ofthe  nanotube  Hence  chiral  nanotubes  require  very  large  fundamental  domains 
and  the  associated  large  computational  cost.  The  use  ofthe  Objective  Structures  framewoik  allows  us  to  consider  very 
small  fundamental  domains,  regardless  of  chirality  with  minimal  computational  effort.  Further,  our  efforts  area  step 
toward  implementing  a  framework  [1,2]  to  understand  the  equilibrium  structures  of  carbon  nanotubes  under  applied 
loads. 
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Mathematical  Foundations  and  Algorithms  for  the 
Quasicontinuum  Method 
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We  will  give  an  overview  ofthe  current  state  of  numerical  analysis  and  algorithms  for  the  quasicontinuum  method  We 
will  present  results  for  atomistic-to-continuum  couplingerror,  rates  ofcorrvergence  for  iterative  solution  methods,  a  pos 
teriori  error  estimation,  and  adaptive  algorithms  for  atomistic-continuum  modeling  and  continuuum  coarsening. 
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Thermoplastic  Elastomers:  Multiscale  Modeling  and 
Microstructure  Evolution 
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■Thermoplastic  elastomers  (TPEs)  are  block  copolymers  made  up  of  "hard’*  (glassy)  and  "soft"  (rubbery)  domains  that 
self  organize  on  a  length  scale  of  a  few  tens  ofnano  meters.  Under  typical  processing  conditions,  TPEs  also  develop  a 
"granular"  structure  at  the  micron  level,  which  is  similar  to  that  of  metal  polycrystals.  Therefor^  TPEs  can  be  described 
as  materials  with  (continuum)  heterogeneities  at  two  different  length  scales.  In  this  talk,  we  will  develop  constitu¬ 
tive  models  for  TPEs  with  lamellar  morphology,  where  the  grains  a  re  made  up  of  the  same,  perfect,  lamellar  structure 
(single  crystal)  with  randomly  varying  lamination  directions  (crystal  orientations).  In  particular,  based  on  experimental 
evidence  we  consider  two  types  ot  such  materials:  "oriented"  and  "unoriented"  samples.  The  oriented  TPEs  are  highly 
ordered,  near  single- crystal  systems  with  slightly  varying  gram  orientations,  while  the  un  oriented  samples  are  initially 
isotropic  dueto  an  initially  random  distribution  ofo  dentations  in  the  sample.  For  oriented  TPEs,  it  is  found  that  under 
certain  loading  conditions;  namely,  under  those  with  sufficiently  large  compressive  deform  at  ions  applied  in  the  direction 
otthelayers  within  the  m d ivi d u al  grain s?t he  overall  behavior  ofthe  material  becomes  macroscop  ically  unstable  (i.e.,  it 
loses  strong  el  lipticity).  The  unoriented  samples  are  also  susceptible  to  loss  of  strong  ell  lpticity,  especially  for  large  con 
trast  in  mechanical  properties  ofthe  hard  and  soft  domains.  Finally  these  list  abilities  will  be  related  to  the  evolution  of 
the  underlying  microstructure,  which  can  be  tracked  experimentally  (via  smallX  ray  scattering  measurements)  in  these 
systems. 
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High  Pressure  Mechanochemistry:  Multiscale  'theory 
Experiment,  and  Search  for  New  Materials 
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High  pressure  mechanochemistry  studies  the  effect  of  large  plastic  shear  on  high  pressure  phase  transformations  and 
chemical  reactions.  We  dev  eloped  rotational  diamond  anvil  cell  (RDAC)  technique  to  study  in  situ  strain  induced 
phase  transformations  in  various  material  systems  [  1  ],  in  eluding  BN,  Si,  SiC  and  iron  Synchrotron  X-ray  diffraction  and 
Raman  spectroscopy  is  used.  A  three-scale  continuum  thermodynamic  theory  and  closed  form  solutions  are  developed 
to  describe  observed  mechanochemical  phenomena  [2,3].  Finite  element  simulations  are  performed  as  well  [4].  At 
the  nanoscale,  a  model  for  strain-induced  nudeation  at  the  tip  of  a  dislocation  pile-up  is  suggested  and  studied.  At  the 
microscale,  a  simple  st  ra  in  *  controlled  kinetic  equation  for  the  strain -induced  phase  transformat  ions  is  thermodynami 
cally  derived.  It  depends  (among  other  parameters)  on  the  ratio  oft  he  yield  strengths  of  phases.  A  macroscale  model 
for  plastic  flow  and  strain -induced  phase  transformations  in  RDAC  is  developed  and  finite  element  solutions  are  found. 
These  models  explain  why  and  how  the  superposition  of  plastic  shear  on  high  pressure  leads  to:  a)  a  significant  (by  a 
factor  of  310)  reduction  ofthe  phase  transform  at  ions  pressure,  b)  reduction  (upto  zero)  of  pres  sure  hysteresis,  c)  the 
appearance  of  new  phases,  especially  strong  phases,  which  werenot  obtained  without  shear,  d)  pressure  growth  during 
phase  transformations  (despite  the  volume  decrease),  ande)  strain  controlled  (rather  than  time-controlled)  kinetics. 
Some  methods  of  characterization  and  controlling  the  phase  transformat  ions  are  suggested  and  the  unique  potential  of 
plastic  straining  to  produce  high-strength  metastable  phases  is  predicted.  Some  new  phenomena  are  predicted  theoreti 
cally  and  confirmed  experimentally  Thus,  phase  transformations  induced  by  rotational  plastic  instability  was  revealed 
and  used  to  reduce  pressure  for  irreversible  phase  transformations  from  hexagonal  to  superhard  cubic  BN  from  55  GPa 
to  5.6  GPa  [5].  Transformation  induced  plasticity  was  rwealed  [1]  which  serves  as  a  positive  mechanochemical  feedback 
for  intensification  of  phase  transformations. 
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A  Quasicontinuum  for  Multilattice  Phase  Transforming 
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The  quasicontinuum  (QC)  method  [1-2]  is  applied  to  materials  possessing  a  multilattice  crystal  structure.  Cauchy-Bom 
(CB)  kinematics,  which  accounts  for  the  shifts  of  the  crystal  basis,  is  used  in  continuum  regions  to  relate  atomic  motions 
to  continuum  deformation  gradients.  To  avoid  failures  otthe  CB  kinematics,  QC  is  augmented  with  a  phonon  stability 
analysis  that  detects  lattice  period  extensions  and  identifies  the  minimum  required  periodic  cell  size  [3].  This  augmented 
approach  is  referred  to  as  “Cascading  Cauchy-Bom  kinematics”.  Applications  to  one  and  two-dimensional  test  problems 
that  highlight  the  salient  features  will  be  presented  In  particular,  an  effective  interaction  potential  material  model  fora 
prototypical  shape  memory  alloy  is  used  to  simulate  a  complete  shape  mem ory  cycle.  The  results  capture  both  tempera 
ture-  and  stress- induced  martensitic  transformations. 
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About  Stress-Induced  Phase  Transformations  in  Shape- 
Memory  Polycrystals 
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Above  the  phase  transform  at  ion  temperature  shape- memory  alloys  show  a  superelastic  behavior  upon  an  applied  stress. 
This  is  due  to  a  phase  transformation  induced  by  the  applied  stress.  It  is  interesting  to  know  the  onset  of  stress -induced 
transformation  from  a  general  perspective  as  well  as  with  respect  to  applications.  Our  aim  is  to  determine  and  study  this 
onset  mathematically  within  the  framework  of  energy  minimization  and  homogenization  theory  [1]. 

We  study  a  model  case  ofscalar  materials,  give  appropriate  definitions  ofthe  transformation  yield  set  in  polycrystalline 
materials  and  show  that  the  Sachs  bound  is  a  sharp  bound  on  the  yield  set.  That  is,  the  Sachs  bound  turns  out  to  charac- 
texizetheonset  of  phase  transformation  in  this  case. 

Furthermore  we  derive  bounds  on  stress-  strain  curves  ofpolycryst  all  me  scalar  shape-memory  alloys;  these  bounds  are 
based  on  an  approach  by  Milton  and  Serkov  [2],  which  was  developed  to  bound  currents  in  nonlinear  conduct  ion  com¬ 
posites.  To  give  specific  examples,  we  consider  four  variant  materials  with  special  textures. 
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Combinatorial  Approach  to  Materials  Discovery 
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I  will  describethe  applications  of  the  combinatorial  approach  forrapidly  optimizing  and  fine  tuning  physical  properties 
of  multifunctional  materials.  In  particular,  we  have  applied  this  technique  to  search  forther  monastic  shape  memory 
alloys  with  minimal  hysteresis  widths.  Synchrotron  microdifFradion  was  performed  on  ternary  composition  spreads  of 
Ni  Ti  Cu  in  order  to  track  the  continuous  change  in  lattice  parameters  across  a  large  region  of  the  phase  diagram  mapped 
on  a  single  Si  wafer.  A  dear  relation  between  the  transformation  matrix  (determined  by  the  lattice  parameters)  and  the 
thermal  hysteresis  of  the  shape  memory  alloy's  were  observed  across  a  large  compositional  region  which  verified  the  non 
linear  theory  of  martmsite 
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Breaching  the  Work  Output  Limitation  of  Ferromagnetic 
Shape  Memory  Alloys 
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In  recent  years,  a  special  class  of  ad  ive  materials  known  as  ferromagnetic  shape  memory  (FSM)  alloys  has  received  much 
attention  due  its  large  actuation  strains  and  fast  response  Despite  ofthe  large  strains,  the  work  output  of  FSM  alloys  are 
limited  by  their  relativelylow  blocking  stress  of  2-6  MPa,  above  which  magnetic  field  induced  strain  is  not  possible  This 
level  ofstress  impedes  the  utilization  of  these  alloys  in  pradical  applications  where  high  work  out  put  by  small  compo 
nents  is  required.  To  date,  attemptsto  enhance  the  work  output  limitation  are  mainly  focused  on  changingthe  concentra 
tions  ofthe  off  stoichiometric  compounds  oftheFSM  alloys,  as  they  set  the  level  ofthe  magnetocrystallineanisotropy 
energy  constant.  In  this  paper  we  suggest  that  the  work  output  limitation  can  be  breached  by  reducing  the  size  ofthe 
actuator,  and  present  an  indication  of  a  possible  fundamaital  relationship  between  size  and  work  output  capabilities. 

Overall,  thebehavior  ofFSM  alloys  involves  motion  oftwin  boundaries  and  is  significantly  influenced  by  its  microstruc¬ 
ture  Based  on  a  theoretical  model,  we  have  shown  that  down-scale  specimens  have  finer  twin  boundary  microstructure 
that  consequently  may  increase  the  blocking  stress.  In  light  of  this,  a  novel  experimental  method  was  realized  to  establish 
this  conjecture  and  to  provide  comp  rehensive  information  on  thebehavior  of  micro  scale  aduators.  A  series  oftests 
demonstrated  no  aduation  strain  redudionupto  extraordinary  loads  of  lOMPa,  which  is  we  believe  to  be  related  to  the 
size  reduction  oftheFSM  alio  yin  our  experiments. 
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Temperature-Induced  Martensitic  Transformations  in  Perfect 
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Solid-to -solid  martensitic  phase  transformations  are  technologically  important  phaiomena  that  result  in  unique 
macroscopic  material  properties  such  as  the  shape  memory  effect  ferromagnetism,  and  ferroelectric  behavior.  In  shape 
memory  alloys,  such  as  CuAlNi  and  NiTi,  the  marten  sit  ^transformation  can  result  from  a  change  in  temperature  or  the 
application  of  stress.  In  fact  both  temperature  induced  and  stress  induced  transformations  are  essential  for  the  existence 
of  shape  memory  behavior. 

An  Effective  Interaction  Potential  model  for  bi-atomic  shape  memory  alloys,  based  on  a  set  of  temperature- dependent 
atomic  pair-potentials,  is  presented.  The  equilibrium  solutions  of  the  govemingnonlinear  equations  are  found,  as  func¬ 
tions  oftemperatureand  applied  stress,  using  symmetry  arguments  and  Branch  Following  and  Bifurcation  techniques. 

To  check  if  a  given  equilibrium  path  is  observable,  its  stability  against  perturbations  ofarbitrary  (with  respect  to  inter 
atomic  distance)  wavelengths  is  investigated.  This  requires  continuum  level  energy  calculations  as  well  as  a  lattice  level 
phonon  spectra  analysis.  Our  work  predicts  the  existenceofa  hysteretic  two -step  temperature  induced  proper  marten- 
s it ic  transformation  from  the  high-temperature  B2  cubic  austenite  phas^  to  an  intermediate  alpha  IrV  orthorhombic 
phase,  to  a  final  B 19  orthorhombic  martensite  phase  Stress-induced  transfer  mat  ion  to  the  alpha -IrV  phase  is  predicted, 
at  high  temperatures,  and  characteristicproperties  such  as  tension -compression  asymmetry  are  captured.  Additionally 
the  transformation  stress  is  found  to  increase  with  mcreasingtemperature  in  agreement  with  experiment.  The  existence 
of  both  temperature-  and  stress -induced  transform  at  ions  indicates  the  possibility  for  shape  memory  behavior.  Finally,  the 
predicted  t  ransformat  ion  parameters  show  good  correspondence  with  experimental  values  for  the  shape  memory  alloys 
CuAl  Ni  an  d  AuCd . 
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Nanoscale  Piezoelectricity 


Pradeep  Sharma 

Department  of  Mechanical  Engineering,  University  of  Houston 

Crystalline  piezoelectricdielectrics  electrically  polarize  upon  application  of  uniform  mechanical  strain.  Inhomogeneous 
strain,  however,  locally  breaks  inversion  symmetry  and  can  potentially  polarize  wen  non  piezoelectric  (centrosym- 
metric)  dielectrics.  Fl exo elect ricty — the  couplingof  strain  gradient  to  polarization —  is  expected  to  show  a  strong 
size- dependency  due  to  the  scalingofstain  gradients  with  structural  feature  size.  In  this  study  using  a  combination  of 
atomistic  and  theoretical  approaches,  we  investigate  the  "effective"  size- depen  dent  piezoelectric  and  elastic  behavior  of 
inhomogeneously  strained  non-piezoelectric  and  piezoelectricnanostructures.  We  aigu$,  through  compjtational  exam 
pies,  the  tantalizing  possibility  of  creating  "apparently  piezoelectric*  nano -composites  without  piezoelect  ric  constituents. 
Further,  to  obtain  analytical  results  and  tease  out  the  novel  physical  insights,  we  analyze  a  paradigmatic  nanoscale  cant  1 
lever  beam.  We  find  that  in  materials  that  are  intrinsically  piezoelect  rig  thefl  exo  electricity  and  piezoelectricity  effects  do 
not  add  lkiearly  and  exhibit  a  nonlinear  interaction.  The  latter  leads  to  a  strong  size- dependent  enhancement  oftheappar 
ent  piezoelectric  coefficient  resulting  in,  for  example  a  "giant"  500%  enhancement  over  bulk  properties  in  BaTiO^  for  a 
beam  thickness  of  5  nm.  Correspondingly  for  non  piezoelectric  materials  also,  the  enhancement  is  non  trivial  (eg.  80  % 
forSnmsize  in  pa  ra  elect  ric  BaTiO^  phase).  Flexo  electricity  also  modifies  the  apparent  elastic  modulus  of  nano  structures, 
exhibiting  an  asymptotic  scaling  of  1/h2  where  "h"  is  the  characteristic  feature  size.  Our  major  predictions  are  verified  by 
quantum  mechanically  derived  force  field  based  molecular  dynamics  fortwo  phases  (cubic  and  tetragonal)  ofBaTiO  . 
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Metallic  foams  and  honeycombs,  with  their  light -weight,  high  specific  stiffness,  and  well  developed  energy  absorption 
characteristics,  are  of  obvious  utility  in  engineering  applications.  However,  these  structures,  often  made  of  aluminum, 
suffer  permanent  deformation  after  crushing.  Cellular  structures  made  from  shape- memory  alloys  (SMAs)  are  particu- 
lady  intriguing  for  theirpotential  to  deliver  shape  memory  and/or  superelasticity  in  a  light-weight  material.  Realization 
ofopai-celled  Nitinol  has  recently  become  possiblevia  a  (newly  discovered  by  Profs.  D.  Grummon  at  Michigan  State 
Univ  and  J.  Shaw  at  Univ.  of  Michigan)  transient  liquid  reactive  brazing  system  for  creating  robust  metallurgical  Nitinol- 
Nitinol  bonds.  With  this  technique,  prototype  sparse  cellular  honeycomb  structures  have  been  made  and  tested,  showing 
up  to  50%  repeatedly  recoverable  strains. 

Of  particular  interest  here,  in  addition  to  modeling  experiments  on  SMAs,  is  an  in  depth  investigation  oftheisother 
mal,  hysteretic  response  of  these  honeycomb  materials  and  thestabilityoftheir  deformation  patterns  at  finite  strains.  A 
detailed  investigation  of  the  influence  ofthe  material  constitutive  law  on  the  macroscopic  stress  strain  behavior  of  these 
honeycombs  is  presented.  It  is  found  that  instabilities — in  the  form  of  departure  from  the  principal  solution  where  all 
cells  deform  in  the  same  pattern — appear  with  the  increase  of  macroscopic  strains  as  the  cell-wall  alloy  transforms  from 
austenite  to  martensite  and  then  eventually  disappear,  upon  further  macroscopic  straining,  as  all  the  cell -wall  alloy  is  all  in 
the  martensit ic phase.  Upon  unloading,  where  the  cell-wall  alloy  material  reverses  this  transformation,  similar  instabilities 
appear  and  disappear.  The  stability  of  honeycomb  structure  is  addressed  i)  byusing  one  cell  Bloch  wave  calculations  for 
the  infinite  perfect  case,  ii)  by  using  large  perfect -geometry  samples  to  account  forspecimen  boundary  conditions  and 
finally  iff)  by  using  large  imperfect -geometry  samples  to  compare  and  contrast  with  experimental  results. 


ABSTRACTS 


12S 


Some  Two-Phase  Problems  of  Nonlinear  Elasticity  in  the 
Presence  of  Small  Interfacial  Energy:  A  Global  Bifurcation 
Approach 
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We  consider  several  different  physical  problems,  having  a  common  mathematical  structure:  a  non-convex  energy  in  the 
lower-order  terms  plus  a  higher-order  regularization  (or  singular  perturbation),  the  latter  characterized  by  a  small  param 
eter  (or  parameters).  We  study  the  equilibria  of  such  systems,  focusingonlocal  minima  of  the  total  energy  We  propose 
a  novel  approach  based  upon  techniques  ofglobal  bifurcation  and  a  priori  bounds — with  the  reciprocal  ofthe  small 
parameter  playing  the  role  of  a  "large*  continuation  parameter  We  indicate  applications  to  shape-memory  solids,  thin 
flexible  sheets  and  lipid  bi-layer  membranes. 
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Ni-Ti-based  Shape  Memory  Alloys 


Dominique  Schryvers* 

Department  of  Physics,  University  of  Antwerp 
mck.schryvers@ua.ac.be 

Remi  Delville 

Department  of  Physics.  University  of  Antwerp 

The  quest  for  new  shape  memory  alloys  better  suited  for  particular  applications  usually  follows  a  path  where  numerous 
extra  constituents  are  added  to  the  basic  alloys  in  smaller  or  larger  amounts.  Though  the  choice  ofelements  and  ratios  is 
normally  based  on  their  known  effects  in  similar  systems,  this  approach  is  essentially  still  trial  and  error  based.  Recently; 
a  new  approach  has  been  suggested  by  R.  James  et  al.  in  which  a  dedicated  search  for  alloys  with  a  low  hysteresis  was 
coupled  with  measurements  of  the  austenite-martensite  transformat  ion  matrix  [1].  Alow  or  zero  hysteresis  is  expected 
when  the  middle  eigenvalue  la  mb  da  2  of  this  tea  ns  for  mat  ion  matrix  equals  1.  At  the  sametime,  this  condition  implies  a 
perfect  fit  at  thehabit  plane  so  that  no  microt  w  inning  in  the  martensite  is  necessary.  Moreover,  this  condition  also  divides 
the  remaining  phase  space  intotwro  parts,  one  with  Type  I  and  II  twins  (lambda!  >  l)  and  one  w  ith  compound  twins 
(lambda2  <  l).  Onesuch  system  is  Ni-Ti-Pd  in  which  Ni  is  replaced  by  Pd. 

The  present  contribution  focuses  on  transmission  electron  microscopy  (TEM)  work  o  f  Ni  Ti  Pd  alloys  with  different 
compositions  with  lambda2  varying  around  1.  The  TEM  micrographs  dearly  confirm  the  existence  ofType  I  and  IT  twins 
at  lambda2  >  1  with  a  decreasing  tw  in  density  when  approaching  lambda2  =  1  In  one  sample,  a  twuiless  habit  plane 
could  be  observed  to  grow  during  an  in  situ  cooling  experiment.  However,  on  some  occasions  also  compound  twins  are 
observed  in  samples  with  lambda!  >  1. 
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Mechanics  of  NiTi  Fatigue  for  Vascular  Stent  Design 


Brian  Berg* 

Senior  Research  Fellow 

R&D  Cardiovascular  Boston  Scientific 

Brian.Berg@bsci.com 

Nitinol  has  become  a  frequently  used  medical  implant  material  despite  its  notorious  complexity  especially  with  respect 
to  fatigue  [  1  ].  This  inter  metallic  compound,  NiTi,  whai  designed  &  manufactured  into  a  femoral  arterial  stent  is  able  to 
spring  out  of  a  int  ra -vascular  catheter,  hold  the  diseased  artery  opai,  and  undergo  millions  of  cycles  of  large  deforma¬ 
tion  due  to  walking,  sitting  standing,  et  cetra  without  failure,  usually  This  talk  will  elucidate  how  mechanics  is  used  to 
measure  and  model  the  use  conditions  in-vivo,  and  the  fatigue  life  of  a  NiTi  stent 
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